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Abstract 

We present the Heisenberg-picture approach to the quantum evolution of 
the scalar fields in an expanding FRW universe which incorporates relatively 
simply the initial quantum conditions such as the vacuum state, the thermal 
equilibrium state, and the coherent state. We calculate the Wightman func- 
tion, two-point function, and correlation function of a massive scalar field. 
We find the quantum evolution of fluctuations of a self-interacting field per- 
turbatively and discuss the renormalization of field equations. 
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I. INTRODUCTION 



Quantum scalar fields such as the inflaton or Higgs fields in an expanding universe are 
characterized by particle production and nonequilihrium state during a rapidly expanding 
period. Path integral and canonical method are two typical methods to describe the quantum 
evolution and oftentimes to incorporate an initial condition. Path-integral formalism requires 
a modification for the nonequilihrium process. The real-time formalism, the imaginary-time 
formalism, and the complex-time formalism are the field theoretic methods to describe such a 
nonequilihrium quantum evolution of the field from an initial thermal equilibrium. The real- 
time formalism has been frequently used to evaluate the finite temperature Green's function 
in the expanding FRW universe The canonical method based on the Hamiltonian 

density operator also has been widely used in which the functional Schrodinger equation 
describes the quantum evolution of the Higgs or scalar field in the expanding universe . 
Functional Schrodinger-picture has been developed to find the wave function, bearing a close 
parallel to the Schrodinger-picture of quantum mechanics. 

The purpose of this paper is to present the Heisenberg-picture approach to describe the 
quantum evolution of the scalar fields in the expanding FRW universe. The Heisenberg- 
picture has the advantage in incorporating an arbitrary initial quantum condition in a par- 
ticularly simple form. For this purpose we reinterpret the quantum operators, the so-called 
generalized invariants 0, for a time-dependent quantum system as constant operators in 
the Heisenberg-picture, and in terms of which we find the quantum field in the Heisenberg- 
picture. Extensive studies have already been done to find and to make use of these invari- 
ants for time- dependent oscillators. A remarkable point is that the exact quantum states 
are the eigenstates of these invariant operators up to time-dependent phase factors in the 
Schrodinger-picture. We applied the generalized invariants to the Schrodinger-picture de- 
scription of the massive scalar field in the expanding FRW universe and Gao et al 
also applied independently the Schrodinger-picture to both the massive scalar field and the 
self- interacting scalar field in Ref. |T^. Recently, the Heisenberg-picture, which has been sel- 
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domly used, however, was employed to find the exact quantum evolution for the same quan- 



tum systems |TT|]. Therefore, it would be worthy to apply the Heisenberg-picture approach to 
the scalar fields in the expanding FRW universe, which can be regarded as an infinite system 
of decoupled time-dependent harmonic oscillators and coupled time- dependent anharmonic 
oscillators for the massive scalar field and the interacting scalar field, respectively. 

II. MASSIVE SCALAR FIELD 

As the simplest model, we consider a massive scalar field in the expanding FRW universe 
with the metric 

ds'^ = -dt^ + a{tfdx^. (1) 



The action is given by 



S = I dtd^x^ 



-AVct^f-im' + inW (2) 



a2 



where ^ = for the minimal coupling and ^ = | for the conformal couphng. We Fourier- 
transform the scalar field 

0(x,t) = ^^/rf^k0k(t)e''^- (3) 
and take the linear combinations 

0ir^w = ^(0k(t)-0-ic(t)). (4) 

Then the x-integration yields 

I dV{x,t) ■.= J2^lit) = I t^'k7rk(t)7r_k(t) (5) 

where a denotes (±,k). Here 0^"^'' and should be treated as independent variables 
corresponding to the cosine and sine modes in the box normalization. Thus one obtains the 
Hamiltonian of mode-decomposed time- dependent harmonic oscillators: 



where 

u^lit) = + ^ + in, (7) 

is the time-dependent frequency squared. Due to the complete mode-decomposition, the 
quantum states is a product of quantum state for each mode: 

i*w)=ni^«w)- (8) 

a 

Each quantum state obeys separately the Schrodinger equation 

lj^\i^a{t))=H^{t)\i;^{t)). (9) 

As pioneered by Lewis and Riesenfeld , the exact quantum states of a time-dependent 
quantum system can be found in the Schrodinger-picture as the eigenstates of the generalized 
invariant operators. We make use of these operators as a method to find the quantum 
evolution of the scalar fields. The invariant operator for the a-th mode harmonic oscillator 
satisfies 

^t-^[J,,^„] = 0. (10) 



There are many operators that satisfy Eq. (|T0|). An invariant operator quadratic in mo- 
mentum and position was found explicitly for a time- dependent harmonic oscillator In 
the Schrodinger-picture, a pair of the invariant operators that are linear in momentum and 
position, and can be interpreted as the time-dependent creation and annihilation operators 
were found P] and were applied to the massive scalar field in the FRW universe How- 
ever, it is our observation that these invariant operators have a simple interpretation in 
the Heisenberg-picture than in the Schrodinger-picture. We note that the time-dependent 
invariant operators in the Schrodinger-picture 

Bi{t) = -i[0a(t)7ra - a^Mt)4>o], 

Bait) = i[Kit)^a - (11) 



where (j)a{t) and 0* (i) are the complex solutions of the classical field equation 

4>a{t) + ^^k{t) + ^l{t)4>a{t) = 0, (12) 

a(t) 

can be interpreted as the time-independent operators in the Heisenberg-picture 

J^BaAt) = ^B^,H - t[Ba,H. Kh] = 0- (13) 

In order to have the usual commutation relation [Ba,H, bI^,h] — 1; the classical solutions 
satisfy the boundary condition 

za3(0,(t)0;(t)-0,(t)0;(t)) = l. (14) 

We denote these by Al^ and Aa, respectively. Thus, we find the position and momentum 
operators in the Heisenberg-picture 

<t>a,H{t) =0a(t)ia + CW^L 
^a,H{t) = AMt)Aa + ^ami]. (15) 

The Fock space of each mode is constructed from the vacuum state 

i„|0„,to) = 0. (16) 
The vacuum state of the scalar field is the product of the vacuum state of each mode 

|0,io) = niO-'^o). (17) 

a 

To apply this Heisenberg-picture to the massive Higgs field and to incorporate the initial 
quantum conditions, we choose Al^ and Aa that diagonalize the Hamiltonian at an initial 
time to 

Ha{to)=u;a{to){AlAa + ^). (18) 
The expectation value of any operator in the Heisenberg-picture is evaluated as 

(O) = {dH)H. (19) 

Below we consider separately three different types of initial conditions. 



A. Vacuum State 

Assuming that the initial quantum state be the nth number state of the Hamiltonian, 
we find the Wightman function of each mode 

{na,to\k{t)Mt')K,to) = n„(0;(t)0,(t') + +0aWC(O- (20) 

From Eq. (^) we find the Wightman function for the scalar field 

(0(x,t)0(x',t'))^™'^ = (O,to|0(x,t)0(x,O|O,to) 



(2vr) 

The two-point function at an equal time follows 



^— I d'ke^^<--'^Mt)<t>m (21) 



x,t)0(x',t))(-^) = 7^/ d'he^'^-^^-^'^lMtW, (22) 



(27r) 

and so does the correlation function 



(02(x,t))(™^) = ^/rf='k|0k(t)p. (23) 



B. Thermal Equilibrium 

Next we consider a thermal equilibrium as an initial state. The initial thermal equilibrium 
is described by the density matrix 

PaAto) = ^e-^^-^^*"), (24) 

where 
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The expectation value of the operator is given by 
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(2sinh^^^)-\ (25) 



ZaJ 



^ T.Me-^^"^'"^dH\n^) (26) 



It is straightforward to evaluate the Wightman function of each mode 



-^lit)Ut') + -. K-lT-,4>.{tWa{t')- (27) 



The Wightman function of the scalar field is obtained by summing over modes 

+ L . te(^)0k(to + K{t)m))\ (28) 



The first term is the vacuum Green's function and the second term is the thermal cor- 
rection 0. The two-point function and correlation function can be found by taking the 
coincident limits t = t' and x = x', t = t', respectively. 

C. Coherent State 

As the last initial condition, we consider the coherent state described by the density 
operator 

p^^^ = ^_g-/3K(io)iJ,iQ+7c^Q+7S^a+ec«]^ ^29) 
Zci.ll 

where = Yia^ joja + co'q,/2. The density operator can be unitarily transformed into the 
canonical form (|^) by acting the displacement operator 

J- * 

D^ = e "<.(to)^"+-<.(to)^" (30) 
and so the partition function becomes 
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We find the Wightman function of the scalar field 

+Mtmt'))] +J^sJ dMKcit)KAt') (32) 



where 

0,,„=«W2!i±MM (33) 

is the classical field corresponding to the coherent state. 

III. SELF-INTERACTING SCALAR FIELD 

We now turn to the case of a self-interacting scalar field. The Hamiltonian takes the 
form 

We consider the potential of a general power law 

yii<p) = Y.^r. (35) 

and the potential of a particular cosmological interest 

VM = ^<P' + f 0^ + ^<P'. (36) 

Unlike the free massive scalar field, the Hamiltonian cannot be simply decomposed as a sum 
over modes, but it is a coupled system of modes. Thus, the canonical approach to the full 
quantum evolution of the fields is extremely difficult to carry out explicitly. However, we 
may find the quantum evolution of fiuctuations around the classical background field and 
include parts of the quantum back-reaction of these fluctuations. By dividing the scalar fleld 
into the classical background field and the fiuctuations 



(37) 



such that 



(0)=0e(t), (^)=0, 



(38) 



and by denoting vr^^ = a^0c, tt^ = a?ip, we rewrite the Hamiltonian as 
i^(t)=/rf^x[(|| + ^(V0.)VaV(0. 

V2a3 ^ 2^^' 2! 502 ^ 



+ (^ + aV0e-V^ + a3 



—77 f + 



(39) 



a-' • 0(/>c ■ 3! 4! 5(/^ 

The canonical formahsm is difficult to apply to the Hamiltonian (|39D , as a whole, consisted of 
the classical background field and the fiuctuations coupled each other. In this paper we treat 
the background field 0c as classical but evolve the fiuctuations (p quantum mechanically. It 
is, however, still difficult to proceed further because we do not yet know the exact quantum 
states for the sub-Hamiltonian of fiuctuations, which is beyond quadratic. 

As a constructive way to find the quantum evolution explicitly, we truncate the Hamil- 
tonian for the fiuctuations at the quadratic order: 



2! 



(40) 



Then we evolve the full fiuctuations perturbatively. Then, as in the massive case, we decom- 
pose the field by modes to the get the Hamiltonian as a sum of time-dependent harmonic 
oscillators now with the frequency squared 



5V(0c 



+ 



(41) 



The quantum back-reaction of the fiuctuations to the classical background field can be eval- 
uated by taking the quantum expectation values of fiuctuation operators in the Hamiltonian: 



i7(0„t) 



'4! 504 



-(^^)(o) 



(42) 



In the above equation the expectation value of any odd power of the momentum and field 
of quantum fluctuations such as {li^^Tifj,), (V0c ■ V(^), ((^), (<^^), ■ ■ ■, vanishes for the initial 
vacuum state and thermal equilibrium state. In these cases we see that 

is an effective potential. The Hamilton equations for the classical field equal to 

0e + 3-0, + -^VefMc) = 0. (44) 



IV. RENORMALIZATION 

It is to be noted that (<^^)(o), ((y5^)(o), etc., contain infinite contributions and require the 
renormalization of the coupling constants A„. Following Ref. ||^, we do this by introducing 
the counter-terms 5A„ in order to cancel the infinite contributions ((^^)(o), (v9^)(o), and etc.. 
The renormalized field equation for the classical background field is 

0. + 3-0.+ (|: 0, +5 in-2)\ '^^ 

\ (rcn) I X \ 

+ E \_4)N r^(^^)(0) + ■ ■ >c = 0, (45) 

where A^'^'^'^^ are the renormalized coupling constants. From the Appendix we find the ex- 
pectation value of quantum fluctuations with the initial vacuum state 

and with the initial thermal equilibrium state 

where each mode of the fluctuations obeys the classical equation 

^„ + 3-^„+(^ + |:^^C->^0. (48) 
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We do explicitly for the potential V2 by introducing the counter-terms 5m?, 5^ and 6\. 
The renormalized field equation for the classical background field reads 

(X 

+3(A(^^'^) + 5A)((^2)(o))0c = O. (49) 

The classical field equations for the fiuctuations involve only the renormalized coupling 
constants 

+ + f ^ + + e(■"^")7^ + 3A(^^°Vc) = 0. (50) 

The Wightman function of quantum fiuctuations at tree level can be found from the 
quantum evolution of ip using the truncated Hamiltonian (^) as in Sec. II. Higher order 
contributions to the Wightman function can be found perturbatively by solving the classical 
background field equation (^) and by substituting the classical background field into the 
Hamiltonian (^H]), and by repeating the procedure. 



V. CONCLUSION 

In this paper we presented the Heisenberg-picture approach to the quantum evolution 
of scalar fields in an expanding FRW universe. The quantum evolution of scalar fields is 
described by the time-dependent functional Schrodinger equation. As a methodology to find 
the wave functions (quantum state) of the functional Schrodinger equation we made use of 
the quantum operators called generalized invariants, which have been used frequently in 
time-dependent quantum systems. These generalized invariants have been applied mostly in 
the Schrodinger-picture. Recently the Heisenberg-picture has been used to find the quantum 



evolution of a time-dependent quantum system [11]. 

The mode-decomposed Hamiltonian of a massive scalar field in the expanding FRW 
universe is a collection of time-dependent harmonic oscillators. The generalized invariants 
that may be used as the creation and annihilation operators have been introduced H and 
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used to find the quantum evolution of tlie scalar field in the Schrodinger-picture 0. In 
the Schrodinger-picture, however, it is rather difficult to incorporate the initial quantum 
conditions. It is our observation that it is relatively easy to incorporate an arbitrary ini- 
tial quantum condition and that the generalized invariants have a simple interpretation as 



constant operators in the Heisenberg-picture |TT|. We found exactly the quantum evolution 



of the massive scalar field incorporating the initial quantum conditions such as the vacuum 
state, the thermal equilibrium state and the coherent state in terms of the classical solutions 
and calculated the Wightman function, the two-point function, and the correlation function 
incorporating the same initial conditions. We also found perturbatively the quantum evolu- 
tion of a self-interacting scalar field by dividing the scalar field into a classical background 
field and fluctuations and by handling the fluctuations quantum mechanically. We found 
the renormalized field equation for the classical background field equation. 
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APPENDIX A: EXPECTATION VALUES 

In this appendix we evaluate the expectation values of the quantum fluctuations with 
respect to the initial vacuum state and the initial thermal equilibrium state. Recollect the 
position operator of each mode of the scalar field 

0a{t) = ^a{t)A^ + ipl{t)Al (Al) 

We put the operator in the normal-ordered form 
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E 



'2n\ f2n - 2\ f2n - 2{k - 1)\ 1 



X : (^„(t)i, + y,:(t)it) 



2n~2k 



(A2) 



to get the vacuum expectation value 



{2n)\ 
2"n! 



(A3) 



In order to evaluate the expectation value with respect to the initial thermal equilibrium 
state, we make use of the following theorem 



(/(a,a^))(*^^™) = (1 -e-")Tr/(a,a^)e-^"^" 



(0, 0\f[Vl + na + Vfic\ Vl + na^ + v^c]|0, 0) 



(A4) 



where 



n 



(A5) 



e'^ - 1 

and c and are the bosonic operators which commute with d and d^, and |0,0) is the 
vacuum state for d and c. Use the theorem to derive the expectation value 



2n\ (therm) 



(0, o| [VTTfi{Mt)Aa + v*a{t)Ai) + Vfi{Mt)ci + ^l{t)Ca)f''\o, 0) 



2n^ 



k=0 



2k 



X [Vfi{^^{t)ci + ^i{t)c^)] |o,o) 



k=0 



2"ra 
{2n)\ 
2"n! 



1 + n)'=n"-'= 



(^:wv^.wr(i+2n)", 



(A6) 



where (7^ and (7], are the bosonic operators which commute with and Aj^,. 
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